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The orthorhombic perovskite NaOsOa undergoes a continuous metal-insulator transition (MIT), 
accompanied by antiferromagnetic order at Tm = 410 K, suggested to be an example of the Slater 
(itinerant) MIT. We have investigated this system using ab initio methods to probe the origin and 
nature of magnetic ordering and the MIT. The rotation and tilting of OsOe octahedra in the GdFeOa 
structure result in moderate narrowing the band width of the tig manifold, but sufficient to induce 
antiferromagnetic (AFM) order. Within the local density spin approximation (LSDA), the system 
magnetically orders but remains metallic with a deep pseudogap. Considering correlation effects 
by including on-site Coulomb repulsion U, at U c ~ 1.5 — 2 eV a MIT occurs only in the AFM 
state. Effect of spin-orbit coupling (SOC) are small in this half-filled t|g shell, although it leads to 
an increase in the critical value U c necessary to open a gap. Our results indicate that this MIT 
is driven by the magnetic order, induced by a combination of structurally-induced band narrowing 
and moderate Coulomb repulsion. Strong p — d hybridization, and not SOC, substantially reduces 
the magnitude of the net moment to ~1 fiB, one third of the formal moment and a value that is 
consistent with the observed moment 

PACS numbers: 71.20.Be, 71.30.+h, 71.27. +a, 75.50.Ee 



INTRODUCTION 



In condensed matter physics the origins of, and 
phenomena associated with, metal-insulator and 
magnetic transitions remain issues of fundamen- 
tal interest, especially when they do not fit into 
the standard mold. From this viewpoint, the oxo- 
osmatcs that display various formal oxidation states 
are attracting great interest, since they display an 
abundant variety of physical phenomena, including 
superconductivity, qualitative changes due to large 
spin-orbit coupling (SOC), unconventional metal- 
insulator transitions (MIT) (the focus of this paper), 
and unusual magnetic behavior. 

In the heptavalent Os 7+ d 1 systems, the cubic 
double perovskite Ba2NaOsOg is ferromagnetic with 
an unusually small magnetic momenti while main- 
taining a cubic structure, whereas the magnetic 
structures of the trigonal NasOsOs, LisOsOe, and 
the double perovskite Ba2LiOs06 are based on an- 
tiferromagnetic (AFM) orderi 2 ^— In the hexavalent 
d 2 systems, the double perovskite Ba 2 CaOs06 is 
antiferromagnetic^ At the metallic d 2 ' 5+ oxidation 
level, KOS2O6 in the /3-pyrochlore structure shows 
superconductivity at T c =8 Kfi- whereas Ba3LiOs20g 
forms a triple perovskite structured This latter com- 
pound becomes AFM at 13K, but shows significant 
reduction in the effective moment of 3.34/ig from 
the theoretical spin-only value of 4.78/is, suggest- 
ing strong SOC effects. In a pentavalent d 3 system, 
the slightly distorted double perovskite Sr 2 CrOs06 



is a semimetallic ferromagnet in which SOC plays 
an important part in the magnetic behavior, by re- 
ducing the spin moment by 0.27/is and inducing an 
orbital moment of 0.1 1[Ib on the Os ion, decreasing 
the net moment by nearly a factor of two £2 Pen- 
tavalent Cd20s207, with its magnetic frustration on 
the pyrochlore Os sublattice, has been suggested to 
be an itinerant Slater insulator^ i.e. a MIT in- 
duced by antiferromagnetic ordering rather than by 
strong correlation effects! 11 ' 12 

The orthorhombic perovskite NaOsOa, also with 
pentavalent Os, synthesized by Shi et al, undergoes 
a MIT coincident with antiferromagnetic ordering 
at the rather high temperature of Tjy=410 K with 
no evidence for any structural distortion below 600 
Kf^ The resistivity and the Hall number decreased 
by ^3 orders of magnitude from Tjv to 200 K and 
continued decreasing to 2K, the insulating charac- 
ter becoming more accelerated below 25 K toward a 
strongly insulating ground state. The electronic con- 
tribution to the specific heat measurement vanished 
at low temperature, indicating a fully gapped Fermi 
surface. Due to this uncommon behavior, seemingly 
a continuously increasing gapping of the Fermi sur- 
face leading to decrease and finally vanishing of car- 
riers, this system was argued to be the first example 
of the continuous MIT (so-called the Slater insulator 
with continuous opening of the AFM gap). 

Recently, Calder et al. confirmed the continuous 
transition as a simultaneous change in magnetic and 
transport character^ The effective Curie moment 
of 2.71 hb observed by Shi et al. is reasonably con- 



sistent with a S = 3/2 moment, reduced (by 30%) by 
hybridization with oxygen and perhaps compensated 
somewhat by an orbital moment. This picture is one 
of high spin d 3 Os 5+ above the transition^ Neutron 
power diffraction (NPD) measurements indicate G- 
type antiferromagnetism (G-AFM) with direction of 
moments alternating on neighboring Os ions, but an 
Os ordered moment of only 1 [Ib^ This reduced 
value is harder to reconcile with with a S = 3/2 pic- 
ture. However, Du et al. calculated a similar (even 
somewhat smaller) value from the sort of first princi- 
ples calculations we will describe below The NPD 
measurements show the spin oriented along the o 
axis without canting of the spins ^ which provides 
additional information related to the impact of large 
spin-orbit coupling (SOC), one aspect that we want 
to investigate more thoroughly. 

Du et al. obtained a magnetic ground state of the 
observed G-AFM type, and provided a comparison 
of this state to other types of magnetic order that 
they could obtain in their calculations for NaOs03^£ 
There remain several unresolved questions: (1) what 
is the origin and the nature of the AFM ordering, 
and what other states are energetically nearby? (2) 
why is the ordered moment only 1 /is for pcntavalcnt 
Os d 3 ? is strong SOC a factor? (3) why is this 
MIT continuous instead of the more common first- 
order type of MIT in transition metal oxides? (4) 
does SOC perhaps play a more central role than Du 
et al. concluded? The orientation of the moments 
(magnetocrystallinc anisotropy) may provide useful 
conditions to address this question. 

For another perovskite-based osmatc Ba2NaOs06, 
we observed strong interplay between correlation ef- 
fects, small bandwidth, and strong SOC that are 
necessary to consider together to build an interpre- 
tation of the unusual features of that system^ which 
include that it is an uncommon ferromagnetic Mott 
insulator, it displays a very small ordered moment, 
and it remains cubic when it seems the d 1 configura- 
tion should favor at least a small Jahn- Teller distor- 
tion. In this paper, we address most of the questions 
we raised above, using first principles spin-polarized 
calculations including both correlation effects and 
SOC. We present work complementing and extend- 
ing that of Du et al. that illustrate more of the 
interplay between correlation effects, SOC, and dis- 
tortion driven band narrowing. 



II. STRUCTURE AND CALCULATION 

Shi et al. reported a highly distorted structure in- 
volving relative displacement of Na and O ions, lead- 
ing to a y/2x 2 x y/2 quadrupled gadolinium orthofcr- 
rite type supercell (space group: Pnma, No. 62).— 




FIG. 1: (Color online) Distorted Pnma crystal structure 
of NaOs03, illustrating the y/2 x v2 doubling cell in the 
a-c plane. There is also a doubling along the b axis due 
to tilting. The red (small) and the yellow (large) spheres 
indicate oxygen and Na atoms, respectively. This figure 
was produced by vesta.— 



TABLE I: Atomic positions of atoms of NaOsOs in the 
Pnma (No. 62) structure, as experimentally observed.— 





site symmetry 
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Na 


4c 


0.0328 
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-0.0065 


Os 


46 
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apical 


4c 


0.4834 


1 


0.0808 


planar O 


8d 


0.2881 


0.0394 


0.7112 



This GdFe03 structure type has been described in 
some detail by Pickett and SinghJ^ The lattice pa- 
rameters arc a=5.3842, &=7.5804, and c=5.3282 (in 
units of A), corresponding to the volume per for- 
mula unit of a cube with edge 3.788 A. The atomic 
positions are given in table [I] This orthorhombic 
distortion leads to three somewhat different Os-0 
bond lengths, but all are close to 1.94 A, resulting 
in a near-ideal aspect ratio of ~T.002iiZ Thus the 
local structure is consistent with the ti g orbitals be- 
ing equally occupied (t 3 g , S=|) retaining local cu- 
bic symmetry in the OsOg octahedron. The O-Os- 
O bond angles are 90.7°, 89.3°, or 89.1°, indicating 
nearly ideal OsOe octahedra. However, the rotation 
and tilting of the octahedra are substantial. The 
O-Os-0 axis tilts by 11° relative to the b axis and 
rotates by 9° in the a (or c) axis. 

In our work both SOC and correlation effects 
have been assessed, within LDA+SOC, LDA+U, 
and LDA+U+SOC approaches, implemented in 
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two all-electron full-potential codes fplo^ and 
wien2Kf^ (The latter combination is available only 
in wien2k.) To include correlation effects in the 
half-filled <2g shell, an effective on-site Coulomb re- 
pulsion U e ff — U — J was used for the results we 
present (except where noted), since experience has 
shown results very similar to those obtained from 
separate U and the Hund's exchange integral J in- 
puts when a shell is half-filled. Convergence was 
checked up to a 14x9x14 fc-mesh of Brillouin zone 
sampling. In wien2k, the basis size was determined 
by RmtKmax = 7, quite adequate for these atoms, 
and the APW radii were Na 2.14, Os 2.15, and O 
1.4, in a.u.. 
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III. ELECTRONIC AND MAGNETIC 
STRUCTURE 

A. LSDA electronic structure 

To begin, we compare electronic structures of the 
cubic and distorted phases. Within LSDA, the cubic 
phase has a higher energy by 378 meV/f.u. than the 
distorted structure in the nonmagnetic (NM) state, 
reflecting the strong lowering of energy by the dis- 
tortion alone. The atom-projected densities of states 
(DOSs) in both phases are compared in the top panel 
of Fig. [5] and are consistent with those presented by 
Shi et alM- As discussed in Ru-based perovskites by 
Mazin and Singh^ the OsOe octahedra must be 
understood as oxygen-shared clusters. In the cubic 
phase, the tig-Pn antibonding (pdir*) and bonding 
(pdir) states lie roughly in the -2.9 to 1.3 cV and 
of -8.7 to -4.7 eV regions, respectively. (The latter 
complex is typically neglected in characterizing the 
electronic structure and in particular the oxidation 
level, leaving the t2g-based pdir* half-filled as the 
physical ti g bands.) Nearly pure p v bands extends 
from -5 to -2.5 eV. The large orthorhombic distor- 
tion leads to narrowing bandwidth and enhancing 
the hybridization. The width of the pdir* manifold 
is reduced by 15%, and a hybridization gap of 1 eV 
appears between the pdir* manifold and the e g -p a 
(pda*) antibonding manifold. 

The bottom panel of Fig. [5] compares the total 
DOS N(E) of the NM cubic, NM distorted, and 
AFM distorted phases, near Ep. In the cubic phase, 
N(£J) varies smoothly very near Ep before giving 
way to a large narrow peak centered 0.3 eV above 
Ep. In the NM distorted phase, a sharp peak ap- 
pears at -60 meV, and van Hove singularities ap- 
pear just above Ep reflecting flat portions of bands 
that may encourage magnetic instability. However, 
our fixed spin moment calculations show no evidence 
of any nearby stable or metastable ferromagnetic 
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FIG. 2: (Color online) Top: orbital-projected densities 
of states (DOSs) for Os 5d and O 2p in the cubic and dis- 
torted Prima crystal structures of nonmagnetic NaOsO:;. 
The distorted phase DOS shows more structure around 
the Fermi energy Ef, reflecting more flat bands that 
enhance magnetic instabilities. Bottom: enlarged total 
DOSs of the nonmagnetic cubic phase and the distorted 
(orthorhombic) nonmagnetic and G-AFM phases, near 
Ef, which is set to zero. Notice that magnetic order 
strongly decreases N(Ef). 



stated In addition to sharp peaks on either side 
of Ep, the AFM DOS of the distorted phase shows 
a deep pseudogap just above E^, suggesting that ap- 
plying U to the Os 5d states at this band filling will 
open a gap. 

Within LSDA, the AFM state with an Os spin 
moment of 0.6 iib has a very slightly lower energy 
of 3 meV/f.u. than NM in the experimentally ob- 
served structure. The local moment of oxygens is 
negligible due to their symmetric placement between 
antialigncd Os moments. This degeneracy (within a 
fraction of fcsTjv/Os) rationalizes the existence of a 



3 





FIG. 3: (Color online) AFM band structure in LSDA 
near Ef (zero of energy) , which is crossed by the Os t2 g 
manifold. A deep and narrow pseudogap is centered at 
Ef- The Os e 9 bands lie higher, above the range of this 
figure, while the O 2p bands lie below -2.5 eV. 

magnetic ordering transition in NaOs03. Figure |3] 
presents the enlarged LSDA band structure, show- 
ing only the Os t2 g (i.e., pdir*) manifold with the 
width W — 3.5 eV. The deep pseudogap in the DOS 
of Fig. [2] results from an incipient gap formation 
arising as a combination of the distorted lattice and 
AFM order. The valence (occupied) band complex 
has become nearly disjoint from the conduction (un- 
occupied) bands in nearly the entire zone, with slight 
overlap remaining at the Y and R points. The top 
two occupied bands are quite flat, accounting for the 
sharp peaks at -0.1 and -0.2 eV as shown in the bot- 
tom panel of Fig. [2] 



B. Origin of the antiferromagnetic order 



FIG. 4: (Color online) Spin density plot of Os ions in 
the a — c plane, with isosurface at 0.13 e/A 3 . This is 
obtained in LDA+U at U e ff = 1 eV, where this system 
becomes insulating, and reflects locally cubic symmetry. 
The different color indicates the opposite spin orientation 
in the G-AFM alignment. 
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FIG. 5: (Color online) AFM total and atom-projected 
DOSs, including primarily the Os t^g manifold, within 
LSDA+U+SOC at U^ff =2.2 eV in the distorted struc- 
ture. In the occupied states below Ef, the unpolarized 
O 2p DOS mirrors closely the majority Os t2 g DOS. The 
result is a greatly reduced spin moment on the Os ion 
(see text). 



To probe the nature and origin of the AFM state, 
we varied both the volume and the angles of the 
OsOg octahedral rotation. In the cubic structure 
with the same volume as the experimentally ob- 
served structure, no AFM order could be obtained. 
To check further for magnetic order in the cubic 
phase, the lattice parameter of the cubic phase was 
varied in the range of ±3%. No magnetic order was 
obtained, indicating that the AFM order is not in- 
duced by simple variation of the Os-Os distance. 
Also, as expected, no AFM is induced by Na dis- 
placement. 

Next we varied the two structural angles individ- 
ually. At a critical tilt of 6 c —3°, a hybridization 



gap between pdir* and pda* manifolds appears, ac- 
companied by AFM ordering with an Os moment of 
0.26 fiB- Somewhat unexpectedly, this moment re- 
mains unchanged up to 9=11°, the observed tilt an- 
gle. When rotating the octahedra by <f>, AFM with 
-Mos=0.13 [Ib appears at <j> c =9°, where the overlap 
between the pdir* and pda* manifolds becomes very 
small. In all cases the opening of a hybridization 
pseudogap around the Fermi level accompanies the 
AFM order, establishing a close connection in this 
compound between structure and magnetic order. 
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FIG. 6: (Color online) Direct and indirect gap vs. effec- 
tive U in LDA+U+SOC calculations. Below the critical 
value of Ueff =2.2 eV, the indirect gap becomes negative. 




FIG. 7: (Color online) Enlarged AFM band structure 
of Os t 2g manifold at U 6f} =2.2 eV in LDA+U+SOC. A 
gap opens at U e ff=l eV in the LDA+U, but at t/ e //=2.2 
eV when SOC is included. This response to f7 e // indi- 
cates that correlation effects are necessary to open a gap, 
however the magnitude of the gap remains an order of 
magnitude small than U e ff- 



C. Metal-insulator transition 

Now we focus more closely on correlation effects 
beyond LSDA on the AFM state in the observed 
distorted structure, by applied the on-site Coulomb 
repulsion U e f / to Os ions using the LDA+U method. 
As might be anticipated from the pseudogap at Ep 
at the required band filling, a gap opens for U e ff as 
small as 1 eV. The Os local moment is 1.0 hb, consis- 



TABLE II: Local moments of Os in the distorted AFM 
NaOs03 in LDA+U+SOC calculations. The moments 
of oxygens are at most 0.01 /is/atom. 



Ueff (eV) 


state 


spin 


orbital 


1 


metallic 


0.69 


-0.050 


2 


semimetallic 


1.03 


-0.085 


2.5 


insulating 


1.13 


-0.097 


3.0 


insulating 


1.22 


-0.106 



tent with the experimentally observed valued This 
value is nearly unchanged when SOC is included (see 
below). Considering 5c? 3 Os 5+ ions, the naive expec- 
tation based on a localized 5d orbital is that orbital 
and spin ordering of ^l^g wou ld be required to ob- 
tain this reduced value of moment. In the insulating 
state, the spin density however assumes a local cu- 
bic symmetry dimpled dice shape, shown in Fig. 21 
which confirms equal occupation of the ti g orbitals. 

Thus the moment must be considered as high spin 
^2g with the moment vastly reduced by 2p hybridiza- 
tion and SOC. Although the reduction in moment is 
common in transition metal compounds and espe- 
cially in the heavy ones,— a reduction by a factor of 
three is highly unusual. 

Figure [6] shows change in the direct and indirect 
gaps due to varying U e ff. After the distortion, ap- 
plication of U e ff, an d SOC, an energy gap opens for 
ll^ff =2.2 eV in this system. Above U^ ff , both gaps 



increase linearly with the difference of 0.3 eV inde- 
pendent of Ueff- Below U£ff, a negative indirect 
gap results due to significant lowering of the bottom 
of the conduction band at the Y point, though the 
direct gap is already considerable (the occupied and 
unoccupied bands are disjoint). This feature may be 
related with the proposal of a continuous transition 
in opening an energy gap, as discussed in another 
osmate Cd2 0s207 by Shinaoka et al^- Considering 
a likely J «0.5 eV, U c 2.7 eV - W suggests for 
this multiband system that a Mott transition is not 
the appropriate picture for this continuous MIT.— 

The corresponding band structure at E^f// ^ s ^is- 
played in Fig. [7] The SOC splitting is about 0.2 eV 
at the T-point. Interestingly, the top of the valence 
bands is nearly flat, reflected in the step increase 
in the valence band DOS that is evident in Fig. [5] 
This indicates that this system is considerably two- 
dimensional. In LSDA+U+SOC, the AFM state is 
nearly degenerate with the NM state, as in LSDA. 
This degeneracy is consistent with the finite temper- 
ature transition between the two states. 
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D. Effects of spin-orbit coupling 

As discussed for several decades, the ti g orbitals, 
when split off from the e g orbitals by the crystal 
field, form a representation of an angular momen- 
tum of L = l<i In the spherical Hund's rule limit 
the total angular momentum is zero for a high spin 
d 3 system. Especially in this locally-cubic geome- 
try and half filling, the orbital moment should be 
small. As shown in Table [Til inclusion of SOC leads 
to an orbital moment antialigned with the spin mo- 
ment and with magnitude no larger than 0.1 /j,b- 
Thus in the insulating state the net moment remains 
about 1 fiB, consistent with the experimentally ob- 
served valued and relatively insensitive to U e ft in 
the range studied here (see Table [TTJ) . 

SOC does affect the band structure in the cubic 
phase, breaking the threefold ti g symmetry down to 
a lower j = \ doublet and a higher j = | quartet. 
For very narrow bands, such as occurs in the double 
perovskite E^NaOsOg, the t^g bands can be split 
into separate j subbands. For NaOs03, however, 
the t2 g bands remain almost 3 eV wide after the dis- 
tortion, thus SOC becomes a secondary effect after 
the distortion, but not necessarily a negligible one. 
We find for example that including SOC leads to a 
change in the bands at the bottom of the insulat- 
ing gap, lowering them in particular at the K-point, 
visible in Fig. [7J 

IV. SUMMARY 

Our results have clarified a number of features of 
NaOs03, and the questions posed in Sec. I. First, a 
combination of the distortion of the perovskite struc- 
ture and moderate correlation effects is responsible 
for both antiferromagnetism and the insulating gap 
in the ground state. Our calculations account natu- 
rally for the factor of three reduction in the observed 
moment (l^s) relative to the ideal t| g value, due to 
very strong hybridization with the O 2p orbitals that 



render the conduction states, and the magnetism, 
itinerant in nature. The oxygen ions, being coor- 
dinated with two Os ions with antialigned spins, 
are unpolarized by (approximate) symmetry The 
calculated energetics, with small energy difference 
between the magnetic and non-magnetic states, is 
also consistent with a Slater-type metal-to-insulator 
transition to an AFM state. Moreover, the simplest 
specification of a Mott insulator is one in which the 
insulating state does not depend on magnetic order. 
In NaOsOa all the evidence indicates that magnetic 
order is necessary for insulating behavior. Because 
of the itinerant nature of the moment, a Hciscnberg 
model cannot provide a representation which could 
then be applied to predict the Neel temperature. 

Spin-orbit coupling, often a dominant force in 5d 
metal oxides, has a more modest effect in this com- 
pound. While it effects the band structure enough 
that the effective coupling strength U e / / is increased 
from 1 eV to 2.2 eV, it is quite ineffective in decreas- 
ing the local moment. The strong moment reduction 
(from a localized high-spin i 2g picture of 3/x_b) is due 
to strong hybridization with O 2p orbitals, which is 
itself affected also by the distortion. The minor im- 
portance of SOC is consistent with the finding of 
Calder et al. that xray resonant photoelectron spec- 
tra show only minor effect of SOCJ^ 

Though 3d perovskite oxides have been studied 
extensively for decades, their 5d counterparts are 
relatively new to details inspection of their prop- 
erties and the underlying mechanisms. Our results 
for NaOs03 will help to provide a basis for under- 
standing the distinctions in the physics of 5d versus 
3d oxides. 
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